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Abstract 



a 

'^ The logarithmic-normal (lognormal) distribution is one of the most frequently ob- 

G served distributions in nature and describes a large number of physical, biological 

Q and even sociological phenomena. The origin of this distribution is therefore of 

'~~' broad interest but a general derivation from basic principles is still lacking. Using 

^— I random nucleation and growth to describe crystallization processes we derive the 

!>■ time development of grain size distributions. Our derivation provides, for the first 
time, an analytical expression of the size distribution in the form of a lognormal 

(T^ type distribution. We apply our results to the grain size distribution of solid phase 

^^ crystallized Si-films. 

o 
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X 

The logarithmic-normal (lognormal) size distribution is one of the most fre- 
quently encountered distributions in nature [HE]- It has been shown early on 
that a large number of technical processes and even sociological phenomena 
such as income distributions [2] or the productivity of researchers [3] follow 
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a lognormal behavior. The distribution is also frequently observed as a re- 
sult of various crystallization processes ^ ; especially random nucleation and 
growth (RNG) processes usually result in a lognormal distribution of cluster 
or grain sizes 0. Lognormal size distributions in phase transformations have 
conventionally been attributed to coarsening [6.7J. Our earlier work |8], how- 
ever, pointed out that such distributions might occur due to time dependent 
kinetics of nucleation and growth, without the involvement of coarsening. Lit- 
erature on nucleation describes the initial or early stages of nucleation PITU] 
and there are a number of derivations of nucleation rates for these early stages, 
see eg. Ref. [H] . A general description of the evolution of size distributions is 
e.g. given in Ref. [12] for aerosols and an overview on the development of 
size distributions in crystallization processes can be found in Ref. [13]. Fur- 
thermore, there are numerical approaches describing the size distributions e.g. 
during condensation of vapor [13] . 

We provide for the first time an analytical derivation for the evolution of log- 
normal size distributions. Preconditions of our derivation are: i. The growth 
starts with the formation of nuclei of a given critical size Qc-, which is usually 
small as compared to the typical size g of grains in the full-grown distribution. 
ii. Nuclei are formed on a random basis during simultaneous growth of other 
grains, which have already grown to a larger size g. iii. There is no coarsening 
(incorporation of smaller grains into larger ones.) These conditions apply to a 
large number of processes termed homogeneous random nucleation and growth 
(RNG). For the case of heterogeneous nucleation, the availability of nucleation 
sites has to be modified by a suitable extension of the present description. We 
first introduce a differential equation that describes the underlying physical 
concepts and present solutions for physically relevant cases. We then apply 
our results to the grain size distribution of solid phase crystallized Si-films. 

We start by establishing the differential equation containing the essential in- 
gredients of a RNG-process. The formation of nuclei requires a critical size gc- 
Therefore, we describe the contribution of nucleation to the grain-size distri- 
bution N at time t by 

^ = /W^te-9.). (1) 

where g is the diameter of grains, J(t) is the time dependent nucleation rate 
and 5 is the Dirac delta-function. Once nuclei are formed, the change of grain- 
size distribution N in time is due to the growth of nuclei and is obtained from 
the continuity equation 

with the particle current density j = A^(5',t) v(g, t). We now assume that the 



growth rate v(g, t) is isotropic and independent of g. Using polar or spherical 
coordinates depending on whether the crystallization of a thin film or a bulk 
material is considered, all angle dependent terms vanish and we obtain 



^''~9'~' dg ■ ^^^ 



This expression holds for dimensions d = 1,2 and 3. [M Introducing Eq. pi) 
into the continuity Eq. ^ and combining Eq. ([I]) with Eq.(|2]), one obtains 
the nucleation and growth equation in d dimensions 



Since the increase in diameter g is determined by the time dependent growth 
rate v{t), the grain size for t —>■ oo may simply be derived from dg{t) = v{t)dt 
and results in 

oo 

ginf = 9c+ v{t)dt. (5) 





We define the quantity N{g,t) = g'^ ^N(g,t) in order to transform Eq. (|4 
into the somewhat simpler shape 



In order to solve this partial differential equation, we use the well-known 
Laplace transform and obtain an ordinary differential equation that can be 
solved analytically. Applying the inverse Laplace transform to the solution 
leads to the final result. The calculations involved in this procedure are be- 
yond the scope of this contribution and will be described separately [15] . 
In order to obtain specific solutions of Eq. Q one has to make a reasonable 
choice for the time dependence of nucleation and growth rates I{t) and v{t). As 
RNG proceeds, nucleation is decreasing due to the decreasing fraction of avail- 
able material, while the average growth rate decreases due to impingement of 
neighbouring grains. The time dependent fraction of material X{t) crystallized 



^ Note, that in contrast to previous work [8], the contribution from the divergence is 
written in terms of the diameter g of the grain rather than its radius r. Consequently, 
the growth-rate v{t) is twice the growth rate in Ref. [8]. 



during RNG-processes is frequently described by the Avrami-Mehl- Johnson 
(AMJ) [TT] expression 



Xd{t) 



1 — exp 



t-toV+' 







t-tr 



(7) 



with the critical crystallization time tc, the incubation time to preceding crys- 
tallization, and the dimensionality of the growth process d = 2 for film growth 
and d = 3 for bulk growth. 0(t) is the Heaviside function. For an experimen- 
tal setup to detect the crystallization kinetics of thin films, see e.g. Ref. [T^ . 
The fraction of material not crystallized at time t and thus still available for 
further nucleation and growth is consequently given by 



Yd{t) = 1 - Xdit) = exp 



t-tr 



d+1' 



e 



t-tr 



In the following we will assume to = 0, as the present paper does not ad- 
dress incubation effects. Following the above discussion, crystallization occurs 
through nucleation and growth processes obeying different time dependencies. 
Accordingly, we introduce two critical time constants td for nucleation and 
tcv for growth. We also consider different, though similar, functional forms for 
the nucleation rate 



/„(t) = /oF„^(t) = /oexp 



tcl 



n+1' 



(9) 



and the growth rate 

Vm{t) = VoY^X^) = Wo exp 



m+l' 



(10) 



and assume m ^ n in order to account for the different processes at work. We 
have found the complete analytical solution of Eq. Q with the above rates. 
Combinations of ra = 1, 2, 3 and m = 1, 2 leads to a variety of functional forms 
that will be described elsewhere [15]. Assuming that nucleation is limited by 
the available material for crystallization, we identify n with the dimension of 
the growth process, and therefore use n = d. The choice of m is, on the other 
hand, less obvious. Choosing m = results i) in a lognormal type distribution, 
and ii) in a more accurate description of the experimental results discussed 
below. The correspondence of our theoretical result with experimental data of 
a specific RNG process may therefore allow to give an independent determi- 
nation of the appropriate value of m, and thereby of the time dependence of 
the growth rate v{t) in the model. In general, Vm{t) will depend on the details 



of the process considered. A description of such details (e.g. the effect of hard 

or soft impingement of grains) is, however, beyond the scope of the present 

work. 

Solving Eq. Q with the rates given in Eqs. (QpO) we obtain for the case m = 



N{g,t) 
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(11) 



(12) 



and the maximal grain size given by Eq. dsl) 

9max{t) =9c + tcvVo (l " 6"*/*=") 



(13) 



Eq. (11) is a remarkable result in that it displays a mathematical structure 



close to the lognormal distribution, a fact not previously recognized in the 
literature. This is seen more clearly when considering the final distribution 
obtained once the nucleation and growth processes have been completed (t — >■ 



oo). Indeed, we then obtain from Eq. (11) 
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(14) 



9inf = lim 9max{t) = 9c + tcvVo- 



(15) 



For 9 ^ 9c ^G can further simplify the result to 
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A comparison of the above result with the lognormal distribution 

2" 



A(^) 



2'nag 



exp 



1 / ln(^/M) y 



with the median M and the width a shows that we obtain a size distribution 
that differs from the lognormal distribution in three respects: a) there is an 
additional dimension-dependent prefactor 1/ g'^~^, b) the exponent of the log- 
normal distribution also depends on the dimensionality of the growth process 
and c) the distribution is multiplied by a cut-off function at the maximal grain 
size {ginf at t ^ oo). Note that Eq. (16) contains the ratio {gd qY which is 



directly proportional to the area or volume of grains. Due to the differences 



between expressions (16) and (18), we cannot directly identify the parameter 
M with ginf or td/tcv with a. 

Although differential equations similar to the one introduced here have been 
discussed in the literature in various fields (see, e.g. Ref. [12], [13] and ref- 
erences therein), our solution departs from those found elsewhere in several 
respects. One of them is the presence of the cutoffs at low and high grain 
sizes in Eqs. (11), (14), (16), which directly result from our derivation and 



the proper treatment of the Dirac distribution. Their appearance correctly 
reproduces the physics of the RNG process. We also note that in the case 
n = 1, m = the particular form of the rates /n(t) and Vmit) exactly leads 
to the lognormal distribution generally introduced empirically to fit data. To 
the best of our knowledge this is the first analytical derivation of a lognormal 
type distribution reported so far. 

We now apply our result to experimental data obtained from the solid phase 
crystallization of amorphous Si as an example for an RNG-process. The experi- 
mentally determined size distribution is taken from Ref. [^ . The corresponding 
parameter fo = 1.72 ± 0.34 /im/h is taken from Ref. [8], which also states the 
crystallization parameters to = 1-5 ± 0.5h and tc = 4 ± 0.5h, necessary for 
Eq. (Is]), using d = 2 for amorphous Si films solid phase crystallized at 600°C. 
The critical grain size gc is of the order of a few nanometers. Since a reliable 
experimental value of /q is not available but Iq and gc only appear as a product 
in the size distribution, we treat the constant, grain size-independent prefac- 



tor of Eq. (16) as a fit parameter. Note that the use of experimental values 
of Jo would have to take into account a factor related to the transformation 
to polar or spherical coordinates when introduced in the solutions described 
above 



Figure 1 compares the grain size distribution of a fully crystallized Si film to 
our calculations above. Fig. la shows a histogram of the distribution of an 
ensemble of about 900 grains. Fig. lb depicts the same set of data, however 
the histogram is created using the distribution of log^^Q g of the individual 
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Fig. 1. Grain size distribution (histogram) of a Si-film with a thickness of 1 ^um 



sohd-phase crystahized at 600C for 15h. Fit to Eq. (16) (sohd fine) and to the 



lognormal distribution given in Eq. (18) (dashed fine) for comparison, a) Distri- 



bution represented on a hnear size scale, b) Distribution represented on a double 
logarithmic scale. 



grains. Every bar of the original distribution of Fig. la represents a number of 
N{g)dg grains in the size interval [g,g + dg\. A transformation of a distribution 
function N{g) to a log-distribution scale is written N{g)dg = g N{\ng)dlng. 
This procedure enables a better visualization of the distribution of small and 
less frequently occurring grains and is applied to the size distribution functions 



derived above 13 



Fitting the lognorinal distribution of Eq. (18) (dashed hne) to the exper 



imentally observed distribution gives M = 1 ± 0.04/xni, a = 0.49 ± 0.03. 



For a fit of the experimental data to Eq. (16) (sohd hne), we use d = 3 
and Qinf = 4.9/im, the largest grain size measured in this distribution. The 
best fit results in td/tcv = 1.74 ± 0.02. Using, as stated above, a growth 
rate of vq = 1.72 ± 0.34/im/h one obtains the following reasonable values 
tcv = 2.9 ± 0.6h and td = 5 ± Ih. It should be noted that tcv and Vq can be 
varied within reasonable limits as long as the ratio td/tcv and g^f remain con- 
stant. A fit using d = 2 was also attempted, but the result is of inferior quality. 
We attribute this finding to the fact that a large fraction of grains have a di- 
ameter smaller than the film thickness of 1/im and growth is therefore mainly 
three-dimensional. Parameters obtained from other crystallized samples lead 
to comparable results. A comparison to further experimental results including 
partially crystallized samples will be presented elsewhere [TB]. We stress the 



fact that, contrary to the fit using the lognormal distribution Eq. (18) (dashed 



line on Fig.l), our fit (solid line) is based on a physical model and theoretical 



derivation of Eq. (16) applied to RNG processes. 



In conclusion, we present a generic differential equation for the time-dependant 
grain size distribution resulting from random nucleation and growth (RNG) 
processes. We specify the nucleation and growth rates according to the well 
established Avrami-Mehl- Johnson-Equation, and for the first time achieve an 
analytical derivation of logarithmic normal-type distributions. We therefore 
propose that the frequently observed lognormal size distribution resulting from 
RNG originates from the dynamics of the process as described here. Our com- 
paratively simple approach is well supported by the grain size distribution 
found in fully crystallized Si-films. Our model may therefore significantly ben- 
efit the understanding and tailoring of technically relevant nucleation and 
growth processes. 



^ The same procedure has been used for Fig. 1 of Ref. [8]. Note that this transfor- 
mation changes the shape of the distribution: while the measured grain-size distri- 
bution with a crystallized fraction of 2.8% is indeed non-monotonic, the calculated 
distribution is not. Although both distributions appear non-monotonic on the loga- 
rithmic size scale, the calculated distribution only becomes non-monotonic once the 
crystallized fraction exceeds about 16%. Nevertheless, the validity of the concept 
presented in Ref. [8] is now supported by the generic derivation given here. 
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Note added in proof 

Recently, we became aware of an early, rarely cited paper of Kolmogorov in 
Dokl.Akad.Nauk. SSR 31, 99 (1941). Using a statistical approach, the author 
derived the logarithmic normal distribution for particle sizes resulting from 
grinding. Our forthcoming paper [15] will discuss the relevance of this approach 
and the resulting distribution of grain sizes obtained. 
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